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q | We show that the observed pulsed optical emission of the anomalous X-ray 

pulsar 4U 0142+61 can be accounted for by both the magnetar outer gap models 
and the disk-star dynamo gap models, therefore is not an evidence in favor of one 
of these models as its responsible mechanism. Nevertheless, the estimated high 
energy gamma-ray spectra from these models have different power-low indices, 
and can be tested by future observations of the Gamma-ray Large-Area Space 
Telescope (GLAST). Furthermore, we show by analytical estimations that the 
expectations of a standard disk model is in agreement with the observed unpulsed 
optical and infrared luminosities of the AXP 4U 0142+61. 
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1. Introduction 

Anomalous X-ray pulsars (AXPs) are neutron stars (NSs) which are characterized by 
their persistent X-ray luminosities (L x ~ 10 34 — 10 36 ergs s -1 ) which are well above their 
rotation powers. Observations show no evidence for the existence of massive companions. 
Their periods are clustered in a remarkably narrow range P ~ 6 — 12 s. Some of the AXPs 
belong to SNRs which show that they are young objects (See Mereghetti et al. 2002, for a 
review of AXPs). The X-ray luminosities, periods, and black body temperatures of AXPs 
are similar to those of soft gamma-ray repeaters (SGRs) (see Hurley 2000 for a review). 
Furthermore, recently reported SGR-like bursts from two of the AXPs (Kaspi & Gavriil 
2002, Gavriil et al. 2002) imply that these sources probably have similar origins. 
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Magnetar models (Dunacan & Thompson 1992, Thompson & Duncan 1995) are suc- 
cessful in explaining the fluences and the luminosities of the normal (L ^ 10 42 erg s -1 ) and 
the giant (L > 10 44 erg s -1 ) super-Eddington bursts of SGRs. In the magnetar models, 
surface magnetic fields are > 10 14 G, and huge magnetic energy release through instabilities 
from inside the NS is the source of these bursts. The persistent thermal X-ray emission was 
suggested to be due to the surface heating by magnetic field decay, while the accelerated 
particles in the magnetosphere by the Alfven waves produced by the small scale fractures on 
the NS surface could be the source of the non-thermal emission (Thompson & Duncan 1996). 
Spin-down torque is provided by the magnetic dipole radiation. Despite the success of the 
magnetar model in explaining the burst energetics, it can account for the period clustering 
of SGRs and AXPs only in one set of decay models with special conditions (Colpi, Geppert 
& Page 2001). 

On the other hand, fall-back disk models, which have difficulties in explaining the short 
time scales of the super-Eddington bursts, can account for the period clustering and the 
persistent X-ray emission of these sources (Chatterjee, Hernquist & Narayan 2000; Alpar 
2001; Marsden, Lingenfelter & Rothschild 2001). The enhancement phase of the persistent 
X-ray emission following the giant flare of the SGR 1900+14 can be explained by the post- 
burst release of the inner disk which has been pushed back by the burst energy (Ertan & 
Alpar 2003) (see also Lyubarsky, Eichler & Thompson (2002) for an explanation based on 
crust cooling after being heated by the burst energy). Moreover, it was shown that the 
young NSs with fall-back disks can produce AXPs within the time scales of the ages of their 
associated supernova remnants with reasonable parameters (Chatterjee & Hernquist 2000; 
Ek§i & Alpar 2003a,b). 

Magnetar models have no detailed predictions for the IR and optical emissions seen from 
the AXPs which, in this work, is one of our main concerns based on the disk models. Out of 
five observed AXPs (and two candidates), four sources were reported to have IR counterparts 
(Kaspi et al. 2002, Israel et al. 2002, Israel et al. 2003). One of them (4U 0142+61) was also 
detected in the R-band (Hulleman et al. 2000). The optical/IR luminosities of these AXPs 
lie above the extrapolation of the black body components of their X-ray spectra (Hulleman 
et al. 2001, Wang & Chakrabarty 2002, Israel et al. 2003) indicating different origins for 
the optical/IR emission, and the thermal X-rays likely to be emitted near the NS surface. 
The observed ratio of the total R-band luminosity to the X-ray luminosity is ~ 1.4 x 10~ 4 
for AXP 4U 0142+61 (Hulleman et al. 2000, Hulleman 2001), and was claimed to be too 
low for the origin of the optical light to be an accretion disk. Kern and Martin (2002) 
reported that the optical (in the R-band) emission from the AXP 4U 0142+61 has a pulsed 
component with a pulsed fraction ~ 0.27 which is much higher than the X-ray pulsed 
fraction ~ 5% (1 — 2 keV), which strongly suggests different emission sites for the pulsed 
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X-ray and optical emissions. It was argued that the pulsed optical radiation is an indication 
that AXPs are magnet ars. 

In present work, we show that the pulsed optical emission can be accounted for by both 
the disk-star dynamo gap model (Cheng & Ruderman 1991) and the magnetar outer gap 
model (Cheng & Zhang 2001), and does not exclude one or the other (Sections 2 & 3). In 
Sec. 4, we show that the two models can be observationally tested through their expected 
different high energy gamma-ray spectral characteristics. Furthermore, in Sec. 5, we show 
by analytical estimations that the R and K-band luminosities expected from a standard disk 
model is in agreement with the observations of the AXP 4U 0142+61. 



2. Disk-Star Dynamo Model 



For simplicity, we consider a NS dipole magnetic field aligned with the spin axis. When 
the inner rim of a Keplerian disk rotates faster than the neutron star, a gap (accelerator) is 
formed around the null surface (f2.B=0) separating the charges in the magnetosphere. The 
electromotive force of the disk-star dynamo (see Cheng and Ruderman (1991) for details, 
CR thereafter) is given by 



£{r)= f (n K (r)-n,) B z {r) — 

Jr C 



(1) 



where B z (r) is the vertical component of the magnetic field through the disk, Vl K (r) the local 
Keplerian angular velocity of the disk, the angular frequency of the NS, and r the inner 
radius of the disk. In a stationary gap, the electromotive force is balanced by the potential 
drop V g across the gap . Substituting Qk(?") = (GM/r 3 ) 1 ^ 2 in Eq.l, and neglecting the effects 
of inflowing disk matter, we obtain 
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where i?*, B*, and M* are the radius, surface magnetic field and the mass of the NS re- 
spectively. Whether the gap is closed or open depends on how close the co-rotation radius 
r c = (GM/fi^) 1 / 3 is to the inner disk radius r . If r is less than 1.8 r c the gap cannot close 
(the geometry of the open and closed gaps are given in Figs. 2 & 3 of CR). In the limit 
fi* <C ^k(^o) the potential drop in Eq.(2) can be as high as 



S(r) ~ 4 x 10 14 rffi B. 
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where we use the notation A = A yi x 10 l . This potential drop can only be reached by a 
static gap in which the pair production is negligible. In a more realistic stationary gap, 
the potential drop V g is limited by the flow of electrons and positrons through the gap, 
especially when there is copious amount of X-ray photons traversing the gap. The X-rays 
are up-scattered by the inverse Compton collisions with the gap accelerated electrons and 
positrons. The gamma-rays produced by these inverse Compton scatterings create secondary 
e ± pairs through the collisions with the X-rays. These secondary pairs emit synchrotron 
emission out of the gap, which we test whether to be responsible for the observed pulsed 
optical emission. For such a pair production limited gap, the maximum potential drop can 
be estimated as 

Vg = 6 x 10 10 L X , 36 (^) " B-tf R'l V (4) 

(CR). The energies of the primary pairs accelerated in this potential drop V g are limited by 
their collisions with the X-ray photons. On the other hand, the gap power is limited by the 
power of the inflowing disk matter from r c down to the inner disk radius r , 

P ~ GMM ( — - — )• (5) 



r 



The inner disk radius can be written as 



1/7 



ro * P f = 2.9 x 10* P £x,36 [jfj cm ( 6 ) 

(Ghosh & Lamb 1979) where r A is the Alfven radius and (5 ~ 1. For the AXP 4U 0142+61, 
r c ~ 7 x 10 8 cm, and the X-ray luminosity is L x = 6.5 x 10 35 ergs s _1 for a source distance 
d = 3 kpc. d was reported to be greater than 2.7 kpc (Hulleman et al. 2000). The accretion 
rate is related to the X-ray luminosity through M = L X R*/GM ~ 3.4 x 10 15 g s _1 . Then, 
Eqs. (5) & (6) give a gap power P g ~ 1 x 10 33 ergs s" 1 . The fraction of the gap power that 
is taken out of the gap by the secondary pairs depends on the pair production optical depth 
of the X-rays to primary gamma-rays. The density of the thermal X-ray photons at r ~ r 
can be estimated as 

(n) L 9 X /rnx ~10 16 cm- 3 (7) 

V 1 47r(r ) 2 c(E x ) V ; 

where (E x ) ps 3kT is the typical energies of the X-rays. For 4U 0142+61, 3kT 1.2 keV 
(Israel et al. 1999, Paul et al. 2000, Mereghetti et al 2002, Juett et al. 2002). The optical 
depth for the pair production of the primary gamma-rays with the soft X-rays becomes 



t ~ in, 



a p I ~ (10 16 cnT 3 ) (10~ 25 cm 2 ) (5 x 10 8 cm) ~ 0.5 (8) 



where a p ps 10 25 cm 2 is the cross section for the pair production, and / ~ r is the length 
of the gap along the field lines (I can exceed r considerably due to the azimuthal winding 
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of the magnetic field lines across the gap). Eq.(8) implies that about half of the gap power 
is transferred to the secondary e ± pairs. We investigate whether the synchrotron emission 
from these secondary pairs just out of the gap can account for the observed optical pulsed 
emission from AXP 4U 0142+61. 

Energy distribution of the secondary e ± pairs can be approximated as 

N e ~ 7 " 2 In (^^j (9) 

(Cheng, Ho & Ruderman 1986) where 7 is the Lorentz factor of the secondary electrons 
and positrons such that 7 m i„ < 7 < 7 max - The synchrotron energy spectrum from these 
secondary pairs is given by 

d 1 NMF(x) (10) 

, min 

where F(x) = x K 5 / 3 (y)dy, K 5 / 3 is the modified Bessel function, x = -E 7 /-E syn , and 
E syn ~ (3/2) ^ 2 h(eB /m e c). The observed energy spectrum becomes 

where AQ is the beaming solid angle of the synchrotron emission, and d is the distance 
to the source. The total synchrotron luminosity L syn is found by integrating f(E^) overall 
frequencies. The lower cut-off of the spectrum is determined by the minimum Lorentz factor 
7min of the secondary pairs for which we take 7 m i n = 1. The maximum Lorentz factor of the 
secondary pairs is 7 max ~ MeV/(E x ) ~ few xlO 3 (Cheng, Ho, & Ruderman 1986; Cheng 
& Ding 1994). The exact position of 7 max does not affect the results significantly, we take 
7max = 3 x 10 3 . The synchrotron radiation is emitted at a radial distance r > ro where the 
magnetic field strength reduces to about ~ 3 x 10 4 G for a dipole field with ~ 10 12 G on 
the surface of the NS. The resultant photon spectrum of the synchrotron emission is given 
in Fig. 1. 

Our aim is to compare the model pulsed optical luminosity in the R-band with that of 
the observations. The ratio of the observed optical luminosity in the R-band (~ 1.7 — 2.1 
eV) to the X-ray luminosity was reported to be w 1.4 x 10~ 4 from the unabsorbed fluxes 
(Hulleman et al. 2000) with a pulsed fraction of the optical light h R ~ 0.3 (Kern & Martin 
2002). Then, the observed pulsed R-band luminosity is L^ bs ~ 2.5 x 10 31 (Afi/4vr) ergs 
s _1 for L x ~ 6 x 10 35 ergs s _1 where Af2 is the beaming angle of the synchrotron emission. 
To compare our results with the observations, we calculated the ratio of the synchrotron 
R-band luminosity to the total synchrotron luminosity and find that L^/L syn w 1 x 10~ 2 . 
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Since the pair production optical depth provided by the soft X-rays to the primary gamma- 
rays is around one half, L syn ~ -P gap /2 ~ 5 x 10 32 ergs s" 1 for the AXP 4U 0142+61, which 
means ^ 5 x 10 30 ergs s _1 . For a beaming angle AQ ~ 2.5, the optical pulsed luminosity 
estimated by the disk-star dynamo model is in agreement with the observations. 

3. Magnetar Model 

The possible gamma-ray emission properties from the outer gaps for both AXPs (with 
soft thermal X-ray input) and SGRs (including the effect of the hard X-ray component) has 
been worked in detail by assuming dipole magnetar fields with B* > 10 14 G (Cheng & Zhang 
2001, Zhang & Cheng 2002). In the outer gap models of AXPs, the primary gamma-rays 
are produced by the curvature emission of the accelerated primary electrons and positrons 
inside the gap (Cheng & Zhang 2001), which is different from the gamma-ray production 
mechanism of the disk-star dynamo model. The primary gamma-ray luminosity, that is, the 
power provided by the magnetar outer gap can be written as 

L 7 ^4.0xl0 32 (J-X M^-tY (-) 4 (-^-) ergss- 1 (12) 
7 V - 5 / V 10 G J V 6s / V 15km / 

(Zhang & Cheng 1997) where / is the fraction of the volume of the outer magnetosphere 
occupied by the outer gap which is given by 

where T 8 (9) is the surface temperature of the magnetar. Its dependence on the polar angle 
9 with respect to the magnetic axis is weak and may be ignored. It should be noted that if 
/ > 1 an outer gap cannot be formed (Zhang & Cheng 1997). For the AXP 4U 0142+61, 
the reported period is p = 8.69 s and the surface magnetic field estimated from B = 3.3 x 
10 19 (PP) 1/2 is around ~ 1.4 x 10 14 G (Mereghetti et al. 2002), and the corresponding 
fractional size of the outer gap is / ~ 0.89. Then the outer gap power of 4U 0142+61 
becomes P g ~ 1 x 10 33 erg s _1 , which is roughly the same as the gap power of the disk-star 
dynamo model. Similar to the disk-star dynamo model gap, secondary pairs, produced by 
the collisions between the primary gamma-rays and X-ray photons originating from the NS, 
emit synchrotron radiation out of the gap. Unlike the disk- star dynamo gap, the outer gap of 
an isolated NS extends to near the light cylinder radius r t = c/fi* which is ~ 4 x 10 10 cm for 
4U 0142+61. Replacing r by ~ rj/3, where the magnetic field decreases to about 10 2 G, and 
setting / = ri in Eqs.(7 & 8), we obtain the pair production optical depth r p ~ 0.03 L x ,36- For 
the maximum primary gamma-ray energies ~ few GeV (Cheng & Zhang 2001), the maximum 
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Lorentz factor of the secondary pairs are ~ few x 10 3 . In the calculations, we set 7 m i n = 1 and 
7max = 3 x 10 3 . Following the same numerical calculations as that of the disk-star dynamo 
model, we found from the resultant spectrum (Eq.10) that the ratio Lr/ L syn is about 3 x 10~ 2 . 
The total synchrotron luminosity of the secondary pairs L syn ~ r p P g ~ 3 x 10 31 ergs s _1 , 
and the corresponding pulsed R-band luminosity is Lr ~ 3 x 10~ 2 L syn ~ 9 x 10 29 ergs s -1 . 
With a solid angle AQ ~ 0.5 of the beamed synchrotron emission , magnetar model can also 
account for the observed pulsed optical luminosity of the AXP 4U 0142+61 in the R-band. 
The expected synchrotron photon spectrum is presented in Fig. 1. The weakness of the 
magnetar model for AXP 4U 0142+61 is that if the surface magnetic field of this source is 
~ 1.4 x 10 14 G (Mereghetti et al. 2002), then the fraction / of the outer magnetosphere 
occupied by the outer gap (Eq. 13) exceeds unity for i?* less than about 13.9 km. This would 
mean that there is no outer gap of 4U 0142+41 at all (Zhang & Cheng 1997). On the other 
hand, if the outer gap indeed exists in 4U 0142+61 this implies that the equation of state 
(EOS) of this source cannot be soft EOS (cf. Shapiro & Teukolsky 1983) with magnetar 
dipole fields. 

The spectral cut-off of the disk-star- dynamo model lies at the X-rays, whereas it is at 
about 10 eV for the magnetar model. However, the expected X-ray luminosity from the 
disk-star dynamo model is much less than the observed pulsed X-ray emission which is very 
probably emitted near the NS, and is not likely to be tested by the observations. In the 
next Section, we show that the two models could be distinguished through their high energy 
gamma-ray spectral properties. 

4. High Energy Emission from the Magnetar and the Disk-Star Dynamo 



The estimated high energy (keV - ~ 10 GeV) emission spectra of the magnetar and 
the disk-star dynamo gap models are different, although their total high energy luminosities 
are similar. The outer gaps of magnetar models, which lie close to the light cylinder, are 
optically thin to the inverse Compton scattering of the X-rays by the gap accelerated pairs. 
These pairs emit beamed curvature radiation along the outer gap. The resultant spectra 
have been worked for both AXPs and SGRs (Cheng & Zhang 2001, Zhang & Cheng 2002). 
The energy distribution of the primary e* 1 pairs in the gap is given by 



(Zhang & Cheng 1997) where 7 D w 2.8 x lO 7 / 1 / 2 ^ 4 P-V^/IO km) 3 / 4 . The Lorentz 
factor 7 is given by 7(2;) ss 7 x~ 3 ^ 4 for x min < x < x max where x — s/ri and s is the local 
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radius of the curvature. For the AXP 4U 0142+61, 7 ~ 5 x 10 7 . The value of x min depends 
on the magnetic inclination angle a, and it is ~ 2/3 for a = 45° and ~ 2/5 for a = 60°. a; max 
could be taken ~ 2 (Cheng & Zhang 2001). The primary pairs with Lorentz factors between 
7min and 7max emit curvature radiation with typical energies E 1 (3/2)fry 3 (c/s) < GeV. 
Above the characteristic energies spectrum has an exponential cut-off, while I(y) oc z/ 1 / 3 for 
E < E y { . The detailed model spectra of AXPs for different period, magnetic field and 
inclination angles are given in Cheng & Zhang (2001). 

On the other hand, for the disk-star dynamo model, the energies of the gap accelerated 
e ± pairs are limited by the inverse Compton collisions by the X-rays. Thereby produced 
gamma-ray photons create e ± pairs via collisions with the X-rays. The energy distribution 
of the pairs in the gap can be approximated by Eq. (9). For a power-law energy distribution 
N e (-f) oc 7~ m of the pairs, the resultant intensity spectrum of the up-scattered X-rays can 
be written as I c (v) oc v~ v with p = (m — l)/2. For the pair distribution given by Eq. (9), 
we have m = 2 + [ ln(7max/7) The variation of the term in square brackets is small 
along the large range of the Lorentz factors. For a variation of 7max/7 from ~ 5 to ~ 10 3 , p 
changes from 0.6 to 0.8. 

We see that the expected high energy gamma-ray spectra of the magnetar and the disk- 
star dynamo models have very different spectral power-law indices. The model luminosities 
are high enough such that they can be tested by the future gamma-ray mission GLAST 
observations (Cheng & Zhang 2001). 



5. Unpulsed Optical/IR Emission 

Since the unpulsed component of the optical/IR emission is expected to be isotropic, 
there is no beaming correction, and therefore (L-^/ L x ) ^ servc d ~ 10~ 4 for the unpulsed emis- 
sion of the AXP 4U 0142+61. The reported IR and optical magnitudes are K m 19.4 mag 
and R ~ 25 mag (Hulleman et al. 2000, Hulleman 2001). The estimated interstellar ex- 
tinction in the R-band is = 4.4 mag, then it is expected that ~ 0.5 in the IR band. 
We estimate that the ratio of the observed K-band flux to that of R-band (-^K/^iOobserved 
is ~ 10 _1 . Now, we will compare these observations with the expectations of a standard 
accretion disk model. 

We take M* = i]M in where i] is the fraction of the disk mass flow rate which is accreted 
on to the NS, and assume that the mass flow through the inner disk regions emitting in the 
K and R bands can be represented by a steady state mass flow. The local dissipation rate 
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D(r) of the disk at the radial distance r in steady state is given by 

q CM 

D(r) c -M ^. (15) 

For L x ~ 10 36 ergs s -1 , the accretion rate onto the NS becomes M* ~5x 10 15 g s _1 ~ i]M in . 
We obtain 

D(r) « 1.2 x 10 14 r 9 3 ?7 _1 ergs s _1 cm -2 (16) 

where rg is the radial distance in units of 10 9 cm. The corresponding local effective temper- 
ature of the disk becomes 

T eff (r) = V4 ~ 3.8 x 10 4 r 9 - s/ V 1/4 (17) 

where a is the Stefan-Boltzmann constant. The central wavelength of the R-band is Ar ~ 138 
nm. The effective temperature that gives the maximum of its black body emission in the R- 
band is kT ~ 1 eV. Equating kT cS = 1 eV in Eq. (17) we find r(kT ~ 1 eV)~ 4.8 x 10 9 ?T 1/3 
cm. Then, 

^ ~ ( R * ^ ~ i x 1(T 4 V -V 3 (18) 

L x \2r(kT = leV) J 1 y 1 

which is in agreement with the observed ratio reported by Hulleman et al. (2000) for the 
values of rj near unity. To estimate the IR flux from the disk, we equate kT eS ~ 0.1 eV in 
Eq. (17), and find r{kT ~ 0.1 eV)~ 1 x 10 2 r 9 ~ lO 11 ^ 1 / 3 cm. Eq. (17) gives (L K /L X ) ~ 
lO -5 ?] -2 / 3 , that is, Lk/-^r ~ 10 _1 (independent of rf) which is also consistent with the 
observations of the AXP 4U 0142+61. On the other hand, shorter wavelength optical and 
UV emission which is expected from the innermost regions of the disk, could be effectively 
absorbed due to increasing IS absorption with decreasing wavelength along these energies. 

Our estimates remain consistent with the observations as long as rj is greater than 
about ~ 30%. Ek§i & Alpar (2003a) showed that fall-back disks can produce AXP period 
and period derivatives in time scales consistent with their SNR ages using dipole magnetic 
fields with ~ 10 12 — 10 13 G. They estimated that rj could be as low as 0.01 for some 
AXPs, and pointed out that r\ could be higher if the mass loss from the disk due to propeller 
effect during the life time of the AXPs and SGRs is taken into account. From the disk model 
parameters of AXP 0142+61 (Eksi & Alpar 2003b), we can estimate its present disk mass 
flow rate as ~ 3 x 10 16 g s™ 1 , that is, rj ~ 15% which is consistent with our rough estimate 
within a factor ~ 2 — 3. More detailed numerical calculations of the optical/IR emission 
from the AXPs and SGRs will be presented in a separate work. 

In present disk models of AXPs, IR and optical luminosities are assumed to be domi- 
nated by the X-ray irradiation of the disk by the central X-rays (Perna, Hernquist & Narayan 
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2000; Perna & Hernquist 2000). For the observed IR emission from the AXP 1RXS J1708- 
40, these X-ray irradiated disk models give unreasonably large inner disk radii (> 10 11 cm) 
and restrict the radial extend of the disk *AR ~ 5 x 10 11 — 10 12 cm) (Israel et al. 2003). 
Optical/IR emission from the low-mass X-ray binary disks has long been known to be X-ray 
irradiation dominated, and this could be the motivation to include the X-ray irradiation in 
the AXP disk models. However, X-ray could be inefficient for the AXP disks due to the 
following reason. The observationally supported irradiation of the outer disk of low mass 
X-ray binaries (LMXBs) has to be indirect, very probably due to a central extended hot 
corona (see e.g. Dubus et al. 1999), because of the self screening of the disk. A central, 
extended hot corona can be fed by the thermally unstable matter from the inner disk where 
the thermal instabilities are most effective (Shaviv and Wehrse 1986). However, the inner- 
most disk regions of AXPs (and SGRs) are cut by the dipole magnetic fields of their neutron 
stars, which are orders of magnitude higher than the magnetic fields of the NSs in LMXBs. 
Therefore, the X-ray irradiation that prevails at the LMXB disks, does not necessarily be 
effective for the AXP and SGR disks, though it is not impossible. Therefore, in our cal- 
culations we ignore the X-ray irradiation. Observed IR/optical luminosity of AXPs do not 
restrict the outer disk radius in our explanation. That is, the outer disk expands freely in 
the low viscosity state, and needs not to be passive (no viscosity) as proposed by Menou et 
al. (2001). On the other hand, we expect that the inner disk radius is around the Alfven 
radius (~ 5 x 10 s cm) which is about an order of magnitude less than the estimated radial 
distance at which most of the optical R-band luminosity is emitted (~ 5 x 10 9 cm). Inner 
rim of the disk can be irradiated by the X-rays, but the resultant emission is expected to 
be in the UV band which is easily absorbed by the ISM. It is hard to estimate the detailed 
emission spectrum from the innermost disk due to the uncertainties in the physics of the 
disk-magnetosphere interaction. Part of the inner disk matter outflowing due to the pro- 
peller effect may also modify the emission from the underlying innermost disk where higher 
energy optical photons emitted. 



6. Conclusions and Discussion 

We have shown that : (i) the observed optical pulsed emission from the AXP 4U 0142+61 
can be explained by means of both the magnetar outer gap model and the disk-star dynamo 
gap model, and therefore does not eliminate either of the models, (ii) the two models, which 
do not provide an observational test through their different synchrotron X-ray luminosity 
due to the dominating X-ray luminosity from the NS, can be tested observationally via their 
estimated high energy gamma-ray spectra by the Gamma-ray Large- Area Space Telescope 
(GLAST), (in) the IR and optical emission expected from a standard accretion disk model is 
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in agreement with the reported unpulsed IR and optical luminosity of the AXP 4U 0142+61 
for the ratio of the accretion rate onto the NS to the disk mass flow rate r\ > 30%. 

Our estimation for rj is roughly consistent with that obtained from the fall-back disk 
model of AXP 0142+61 (Ek§i & Alpar 2003b). The IR to X-ray luminosity ratio of the 
AXP IE 2259+58 is similar to that of 4U 0142+61 (Hulleman 2001), and can be explained 
by similar arguments. The reported L R /L X ratio of AXP 1RXS J1708-40 is ~ 2 x 10~ 3 
(Israel et al. 2003) which is about 10 times higher than that of 4U 0142+61, and might 
indicate relatively more efficient mass outflow due to the propeller effect in this system with 
i] ~ 1/30. Detiled numerical disk calculations to estimate the optical/IR luminosities of 
AXPs and SGRs and comparisons with their period and period evolution models will be 
presented in future work. 
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Fig. 1. — Photon flux spectrum of the synchrotron radiation from the magnetar outer gap 
model (dashed line) and the disk-star dynamo model (continuous line). The vertical lines 
show the borders of the R-band (1.7 - 2.1 eV). 



